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Abstract. Theory of trembling motion [Zitterbewegung (ZB)] of charge carriers 
in various narrow-gap materials is reviewed. Nearly free electrons in a periodic 
potential, InSb-type semiconductors, bilayer graphene, monolayer graphene and 
carbon nanotubes are considered. General features of ZB are emphasized. It is 
shown that, when the charge carriers are prepared in the form of Gaussian wave 
packets, the ZB has a transient character with the decay time of femtoseconds in 
graphene and picoseconds in nanotubes. Zitterbewegung of electrons in graphene 
in the presence of an external magnetic field is mentioned. A similarity of ZB in 
semiconductors to that of relativistic electrons in a vacuum is stressed. Possible 
ways of observing the trembling motion in solids are mentioned. 



PACS numbers: 73.22.-f, 73.63.Fg, 78.67.Ch, 03.65.Pm 
1. Introduction 

Zitterbewegung (the trembling motion) was theoretically devised by Schroedinger [TJ 
after Dirac had proposed his equation describing free relativistic electrons in a vacuum. 
Schroedinger showed that, due to a non-commutativity of the quantum velocity with 
the Dirac Hamiltonian, relativistic electrons experience the Zitterbewegung (ZB) even 
in absence of external fields. The frequency of ZB is about uj = 2moc 2 /h and its 
amplitude is about the Compton wavelength A c = h/moc w 3.86 x 10~ 3 A. It was later 
understood that the phenomenon of ZB is due to an interference of electron states with 
positive electron energies (E > m c 2 ) and those with negative energies (E < m c 2 ). 
In other words, the ZB results from the structure of the Dirac Hamiltonian, which 
contains both positive and negative electron energies. It is a purely quantum effect 
as it goes beyond Newton's first law. To our knowledge, the ZB for free electrons has 
never been directly observed. However, in the presence of the Coulomb potential the 
ZB is manifested in appearance of the so called Darwin term. It was pointed out that 
the Zitterbewegung also may occur in non-relativistic two-band systems in solids [2]. 
Since the appearance of papers by Zawadzki [3. and Schliemann et al [4] the subject 
of ZB became popular and it was demonstrated that this phenomenon should occur 
in various situations in solids HI El HI [101 [Til H21 H3] . It is the purpose of the 
present review to outline main features of ZB in narrow gap semiconductors. 

We begin by elementary considerations based on the Schroedinger equation. In 
the absence of external fields the Hamiltonian is H — p 2 /2m and the velocity is 
ij = dH/dpj = pj/m. The time derivative of velocity is easily calculated to give 
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Vj = l/(ih)[iij, H] = 0. This means that Vj(t) = const, which is equivalent to first 
Newton's law stating that in absence of external forces the velocity is constant. 

The situation is different when a Hamiltonian H has a matrix form, like in 
the Dirac equation for relativistic electrons in a vacuum or in case of two or more 
interacting energy bands in solids. Then the quantum velocity Vj = dH/dpj, which is 
also a matrix, does not commute with the Hamiltonian and the quantum acceleration 
dvj/dt does not vanish even in the absence of external fields. Below we consider such 
situations and investigate their consequences for the electron motion. 



2. Theory and results 



We consider first the case of InSb-type narrow-gap semiconductors (NGS), see [3J. 
Their band structure is well described by the model of r 6 (conduction), T 8 (light 
and heavy hole) and Ti (split-off) bands and it represents an 8 x 8 operator matrix. 
Assuming the spin-orbit energy A>£ 5 , neglecting the free-electron terms and taking 
the momentum components p z ^ and p x — p y = 0, one obtains for the conduction 
and light hole bands the Hamiltonian 
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where a 3 and (3 are the well known 4 x 4 Dirac matrices and u = (Sj^mJ) 1 / 2 1x10 s 
cm/s is the maximum velocity. Hamiltonian ([1]) has the form appearing in the Dirac 
equation. The electron velocity is z = (l/ih)[z, H] = ua 3 . To determine aa(t) one 
calculates the commutator of a.3 with H and integrates the result with respect to time. 
This gives 2(t), and z(t) is calculated integrating again. The result is 
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where Aq = d(0) — up z /H. The first two terms represent the classical motion. The 
third term describes time dependent oscillations with the frequency ujz ~ 8 g /h. Since 
~ 1, the amplitude of oscillations is hu/H rs h/(m,Qu) = \z- Here 

Az = -Ar (3) 

TTIqU 

is an important quantity analogous to the Compton wavelength A c = h/(m^c) for 
relativistic electrons in a vacuum. The oscillations analogous to those described in 
^ are called Zitterbewegung. The quantity X z can be measured directly. The 
energy for electrons in NGS can be written in a 'semi- relativistic' form (see |14j ) 



E 



k 2 ) 1 / 2 . For k 2 > this formula describes the conduction and light- 
hole bands. For imaginary values of k there is k 2 < and the above formula describes 
the dispersion in the energy gap. This region is classically forbidden but it can become 
accessible through quantum tunnelling. Figure [T] shows the data for the dispersion in 
the gap of InAs, obtained from tunnelling experiments. The fit gives Xz ~ 41. 5A and 
u s» 1.33 x 10 s cm/s, in good agreement with the estimation for InAs (mj ~ 0.024rno)- 
Similar data for GaAs give Xz between 10 and 13 A Q~5], again in good agreement 
with the theoretical predictions (mj w 0.066m ). 

To demonstrate universality of the two-band situation we consider the well known 
case of nearly free electrons in a solid in which the periodic lattice potential V(r) is 
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Figure 1. Energy- wave vector dependence in the forbidden gap of InAs. Various 
symbols show the experimental data of Parker and Mead (Ref. |16|). The solid line 
is a theoretical fit. The determined parameters are \ z = 41.5Aand« = 1.33 x10 s 
cm/s (after Ref. [3]). 



treated as a perturbation (see [9]). Near the Brillouin zone boundary the Hamiltonian 
has, to a good approximation, a 2 x 2 form 

*-(*!? *')■ (4) 

where V* = V- q are the Fourier coefhcients in the expansion of V(r), and = 
h 2 k z /2mo is the free electron energy. The 2x2 quantum velocity v z can now be 
calculated and the acceleration v z is computed in the standard way. Finally, one 
calculates the displacement matrix £y. 

Until now we treated the electrons as plane waves. However, Lock |17j in his 
important paper observed that such a wave is not localized and it seems to be of a 
limited practicality to speak of rapid oscillations in the average position of a wave of 
infinite extent. Since the ZB is by its nature not a stationary state but a dynamical 
phenomenon, it is natural to study it with the use of wave packets. These have become 
a practical instrument when femtosecond pulse technology emerged (sec Ref. [18 ). 
Thus, in a more realistic picture the electrons are described by a wave packet 

1 d 1 / 2 r°° f 1 \ f 1 



ip{z) = -j=—yji J exp y~-^d (k z - k zQ ) J exp(ik z z)dk z y Q 

In Figure [2] we show the ZB oscillations of zu{t) averaged over wave packet ([5]). 
The essential result is that, in agreement with Lock's general predictions, the ZB 
oscillations of a wave packet have a transient character, i.e. they disappear with time 
on a femtosecond scale. 

Now we turn to interesting and intensively studied materials: bilayer and 
monolayer graphcnc and carbon nanotubes (CNT), see pJQIJlIIS]. Two-dimensional 




Figure 2. Transient Zitterbewegung oscillations of nearly- free electrons versus 
time, calculated for a very narrow wave packet centered at various k z o values. 
The band parameters correspond to GaAs. 



Hamiltonian for bilayer graphene is well approximated by [20] 

ft - 1 f (Px-iPy) 2 \ ,(.-. 

HB --2^\{p x+l Pyf J' (6) 

The energy spectrum is £ = ±h 2 k 2 /2m*, i.e. there is no energy gap between the 
conduction and valence bands. The position operator in the Heisenberg picture is a 
2x2 matrix x(t) = exp(ii?Bi/7i):r exp(— lifei/ft). We calculate 
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where k 2 = k 2 + k 2 . The third term represents the Zitterbewegung with the frequency 
hujz = 2h 2 k 2 /2m* , corresponding to the energy difference between the upper and 
lower energy branches for a given value of k. In order to have the ZB in the x 
direction one needs a non-vanishing momentum k y in the y direction. If the electron 
is represented by a two-dimensional Gaussian wave packet similar to that given in 
Equation ([5]) , the integrals of interest can be calculated analytically. In Figure [3] we 
show observable physical quantities related to ZB. It can be seen that they have a 
transient character. 

Looking for the reason of transient character of ZB the electron wave packet is 
decomposed into sub-packets corresponding to positive and negative electron energies. 
It turns out that the 'positive' and 'negative' sub-packets move in the opposite 
directions with the same velocity v = hko y t/2m*, where hko y is the initial value 
of momentum. The relative velocity is v rel = hko y t/m*. Each of these packets has 
the initial width d and it spreads (slowly) in time. After the time T^ 1 = d/v rel the 
distance between the two packets equals d, so the sub-packets cannot interfere and the 
ZB amplitude diminishes. This reasoning gives the decay constant Tz = hko y /m*d, 
in agreement with the analytical results [TO] • Thus, the transient character of the ZB 
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Figure 3. Zitterbewegung of a charge carrier in bilayer graphene versus time, 
calculated for a Gaussian wave packet width d = 300A and ko y = 3.5 X 10 8 m — 1 : 
a) displacement, b) electric current. The decay time is about T~ = 40 fs (after 
Ref. [10]). 



oscillations in a collisionless sample is due to an increasing spacial separation of the 
sub-packets corresponding to the positive and negative energy states. 

Now, we turn to monolayer graphene. The two-dimensional band Hamiltonian 
describing its band structure is [ST] 

Hm = u( . ° (8) 

V Px + Wy J 

where u sa 1 x 10 8 cm/s. The resulting energy dispersion is linear in momentum: 
E = ±uhk, where k — \Jk^. + ky. The quantum velocity in the Schroedinger picture is 

rij = 8Hm I 'dpj, it does not commute with the Hamiltonian. In the Heisenberg picture 
we have v(t) = exp(iHMt/h)v ex.p(—iHjMt/K). Using Equation {8J we calculate 

vi n) =u^sm(2ukt). (9) 

The above equation describes the trembling motion with the frequency uz — 2uk, 
determined by the energy difference between the upper and lower energy branches 
for a given value of k. As before, the ZB in the direction x occurs only if there is a 
non- vanishing momentum hk y . The results for the current j x = ev x are plotted in 
Figure [4] for koy = 1.2 x 10 9 m _1 and different packet widths d. It is seen that the 
ZB frequency does not depend on d and is nearly equal to u>z, as given above for the 
plane wave. On the other hand, the amplitude of ZB does depend on d and we deal 
with decay times of the order of femtoseconds. 

Finally, we consider monolayer graphene sheets rolled into single semiconducting 
carbon nanotubes (CNT). The band Hamiltonian in this case is similar to Equation 
|8| except that, because of the periodic boundary conditions, the momentum p x is 
quantized and takes discrete values hk x = Kk nv , where k nv = (2ir/L)(n — u/S), 
n = 0, ±1,..., v = ±1, and L is the length of circumference of CNT [22]. As a 
result, the free electron motion can occur only in the direction y, parallel to the tube 
axis. The geometry of CNT has two important consequences. First, for v — ±1 there 
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Figure 4. Oscillatory electric current in the x direction caused by the ZB 
in monolayer graphene versus time, calculated for a Gaussian wave packet with 
k 0y = 1.2 X 10 9 m- 1 and various packet widths d (after Ref. [10]). 



always exists a non- vanishing value of the quantized momentum Kk nv . Second, for 
each value of k nv there exists = — k nv resulting in the same subband energy 



The time dependent velocity v y (t) and the displacement y(t) can be calculated for 
the plane electron wave in the usual way and they exhibit the ZB oscillations (see Ref. 
[7]). For small momenta k y the ZB frequency is htoz — E g , where E g = 2%uk nil . The 
ZB amplitude is Xz ~ l/k nv . However, we are again interested in the displacement 
y(t) of a charge carrier represented by a one-dimensional wave packet analogous to 
that described in Equation © , see Ref. fTUj ■ Figure [5] shows the ZB oscillations 
calculated for a Gaussian wave packet of two widths. The decay times are of the order 
of pikoseconds, i.e. much larger that in bilayer and monolayer graphene. The reason 
is that the ZB oscillations occur due to the 'built in' momentum k x arising from the 
tube's topology. In other words, the long decay time is due to the one-dimensionality 
of the system. 

The last subject we consider is the trembling motion of electrons in monolayer 
graphene in the presence of an external magnetic field B \\ z, see [II] . The magnetic 
field is known to cause no interband electron transitions, so the essential features 
of ZB, which results from an interference of positive and negative energy states of 
the system, are expected not to be destroyed. On the other hand, introduction of 
an external field provides an important parameter affecting the ZB behavior. The 
essential feature introduced by the magnetic field is a quantization of the electron 
spectrum E ns = sTiuj^fn, where n = 0, 1, . . . and s = ±1 for the conduction and valence 
bands, respectively. The basic energy is Tito = ^/2hu/L, with 1/L = (eB /K) 1 / 2 . The 
velocity operators can again be calculated in the time-dependent Heisenberg picture. 
These are subsequently averaged over a Gaussian wave packet and they exhibit the 
trembling motion. The presence of a quantizing magnetic field has very important 
effects on ZB: (1) For B ^ the ZB oscillations are permanent, for B = they are 
transient. (2) For B ^ many ZB frequencies appear, for B = only one frequency 
is at work. (3) For B ^ both interband and intraband (cyclotron) frequencies 



£ = ±E, where E = hu^ k 2 nl , + k 2 
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Figure 5. Zitterbewegung of charge carriers in the ground subband of a single 
carbon nanotube of L = 200 A versus time (logarithmic scale), calculated for 
Gaussian wave packets of two different widths d and ko y = 0. After the ZB 
disappears a constant shift remains. The two carriers are described by different 
quantum numbers v (after Ref. |10j 1. 



contribute to ZB, for B = there are no intraband frequencies. (4) Magnetic field 
intensity changes not only the ZB frequencies but the entire character of ZB spectrum. 

3. Discussion 

It follows from the recent theoretical papers that the phenomenon of ZB should 
appear quite often in solids. Whenever one deals with interacting energy bands, 
charge carriers having a non-vanishing momentum should experience the trembling 
motion. Its frequency ujz is about fouz ~ A£, where AS is the energy separation 
between the interacting 'positive' and 'negative' states, and the amplitude Xz goes 
as Xz oc 1/A£. The latter can reach the values of tens of nanometers in narrow- 
gap semiconductors. As mentioned in the Introduction, a similar phenomenon should 
occur for free relativistic electrons in a vacuum, but both its frequency and amplitude 
are much less favorable than for electrons in semiconductors, see Refs. [23] 121] ■ This 
correspondence of the trembling motions illustrates a general analogy between the 
behavior of electrons in narrow-gap semiconductors and that of relativistic electrons in 
a vacuum [21]. The Zitterbewegung was also proposed for other systems [231 [2B], and it 
has some interesting analogies in quantum optics [27] . As to the ways of observing the 
ZB, one can try to detect directly the moving charge at corresponding frequencies using 
the scanning probe microscopy [28[ 129] . One can also try to measure the ac current 
related to the oscillating charge, see Figures [3] [4] and Ref. [H]). Finally, and this is 
probably the most promising way, the oscillating charges should produce an observable 
dipole radiation. An external magnetic field can continuously change the frequency 
of such radiation. It appears that graphene in the presence of a magnetic field 
provides most favorable conditions for an observation of the fascinating phenomenon of 
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Zitterbewegung. Very recently, an acoustic analogue of Zitterbewegung was observed 
in a two-dimensional sonic crystal [30j . 
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